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Abstract —  Polarization  diversity  has  proved  to  be  a  useful  tool 
for  radar  detection,  especially  when  discrimination  by  Doppler 
effect  is  not  possible.  In  this  paper,  we  address  the  problem  of 
improving  the  performance  of  polarimetric  detectors  for  targets 
in  heavy  inhomogeneous  clutter.  First,  we  develop  a  polarimetric 
detection  test  that  is  robust  to  inhomogeneous  clutter.  We  run 
this  polarimetric  test  against  synthetic  and  real  data  to  assess  its 
performance  in  comparison  with  existing  polarimetric  detectors. 
Then,  we  propose  a  polarimetric  waveform-design  algorithm  to 
further  improve  the  target-detection  performance.  A  numerical 
analysis  is  presented  to  demonstrate  the  potential  performance 
improvement  that  can  be  achieved  with  this  algorithm. 


I.  Introduction 

The  detection  of  static  or  slowly  moving  targets  in 
heavy-clutter  environments  is  considered  a  challenging  prob¬ 
lem,  mainly  because  it  is  not  possible  to  discriminate  the 
target  from  the  clutter  using  the  Doppler  effect.  Polarization 
diversity  provides  additional  information  that  enhances  the 
detection  of  targets,  particularly  under  the  conditions  described 
above.  Detection  performance  could  be  further  improved  if  the 
polarization  of  the  transmitted  signal  were  optimally  selected 
to  match  the  target  polarimetric  aspects. 

Earlier  work  in  the  field  of  polarimetric  detection  has 
addressed  the  design  of  detectors  under  the  assumption  that 
training  data  is  available  [l]-[3].  The  performance  of  these 
detectors  can  be  severely  degraded  in  the  presence  of  in¬ 
homogeneous  and  non- stationary  clutter.  To  overcome  this 
problem,  the  application  of  compound- Gaussian  distributions 
for  modeling  the  global  behavior  of  nonhomogeneous  clutter 
has  been  proposed.  However,  the  use  of  non-Gaussian  models 
increases  the  difficulty  of  developing  efficient  detection  algo¬ 
rithms.  For  instance,  the  detection  statistic  derived  from  the 
scheme  presented  in  [4]  has  a  closed-form  expression  only  for 
the  special  case  of  two  polarimetric  channels,  and  the  detector 
proposed  in  [5]  does  not  support  the  constant  false-alarm  rate 
(CFAR)  property  with  respect  to  the  clutter  covariance. 

In  this  paper,  we  first  develop  a  polarimetric  detector 
whose  decision-making  capability  is  based  only  on  the  data 
collected  from  the  range  cell  under  test,  without  resorting  to 
secondary  data  or  prior  knowledge  of  the  target  and  clutter. 
This  feature  makes  it  robust  against  inhomogeneous  and 
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non-stationary  clutter.  Later  in  this  paper,  we  address  the 
problem  of  waveform  design  by  presenting  an  algorithm  for 
selecting  the  waveform  polarization.  The  developed  test  statis¬ 
tic  has  a  closed-form  expression  that  incorporates  information 
about  the  transmitted  signal  polarization.  Hence,  analysis  of 
detection  performance  is  used  to  select  the  next  transmission 
polarization,  which  improves  the  target  detection. 

II.  Polarimetric  Radar  Model 

We  consider  a  mono- static  radar  capable  of  transmitting 
waveforms  with  any  arbitrary  polarization  on  a  pulse-by -pulse 
basis.  In  addition,  we  assume  that  the  radar  system  illuminates 
a  point  target.  The  recorded  data  from  the  range  cell  under  test 
consist  not  only  of  the  target  echoes  but  also  of  the  undesired 
reflections  from  the  target  environment  (clutter). 


A.  Polarimetric  Data 

The  output  of  a  diversely  polarized  array  of  Q  sensors 
receiving  the  echoes  from  the  cell  under  test  can  be  expressed 
as 

y(t)  =  5(5*  +  Sc)  £(t)  +  e(t),  t=l,...,N,  (1) 


where 

•  The  Q  x  1  vector  y(t)  is  the  complex  envelope  of  the 
measurements. 

.  The  Q  x  2  matrix  B  is  the  response  of  the  diversely 
polarized  sensor  array.  If  the  receiver  array  is  a  vector 
sensor  [6],  the  array  response  is  given  by 
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where  p  and  p  are  the  azimuth  and  elevation  angles  of  the 
cell  under  test,  respectively.  For  a  conventional  polarized 
radar  measuring  the  horizontal  and  vertical  components 
of  the  electric  field  and  assuming  these  two  sensors  are 
orthogonal  to  the  direction  that  points  toward  the  cell 
under  test,  the  array  response  matrix  is  B  m  I2. 
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•  The  complex  scattering  matrix  S  represents  the  polariza¬ 
tion  change  of  the  transmitted  signal  upon  its  reflection 
on  the  target  or  clutter: 
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where  the  variables  s n  and  822  are  co-polar  scattering 
coefficients  and  812  and  821  are  cross-polar  coefficients. 
For  the  mono-static  radar  case,  812  =  821.  The  super¬ 
scripts  t  and  c  refer  to  the  target  and  clutter. 

•  The  vector  £(£)  is  the  narrowband  transmitted  polarized 
signal  which  can  be  represented  by 
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where  £1  and  £2  are  the  signal  components  on  the 
polarization  basis  of  the  transmitter,  a  and  (3  are  the 
orientation  and  ellipticity  angles,  respectively,  and  s(t) 
is  the  complex  envelope  of  the  transmitted  signal. 

•  The  vector  e(t)  represents  the  thermal  noise. 

Equation  (1)  can  be  written  as 


y(t)  =  s(t)B£(n  +  x)+e(t),  (5) 


where  the  scattering  coefficient  vectors  of  the  target  and 

clutter  are  y  =  Kh>  4v>shvF  and  x  =  Kh>  4v>  shvF’ 

respectively,  which  have  dimension  P  =  3.  The  polarization 
matrix  £  is 


c  = 
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The  time  samples  can  be  stacked  in  one  vector  of  dimension 
NQ  x  1: 

y  =  (s  <g)  B£)  (y  +  x)  +  e,  (7) 

where  s  =  [s(l s(N)]T  and  ®  is  the  Kronecker  product. 
Piling  together  the  data  corresponding  to  pulses  of  different 
polarization  yields 


y  =  Ay.  +  Ax  +  e, 


(8) 


where  A  is  a  M  x  P  ( M  =  KNQ)  complex  matrix  that 
represents  the  system  response: 

r  s  <g>  sci  l 
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[s  eg.  B^k\ 

and  is  the  polarization  matrix  of  each  diversely  polarized 
pulse  (k  =  1, . . . ,  K). 


B.  Statistical  Model 

We  assume  that  the  target  is  a  small  man-made  object; 
hence,  /x  is  a  deterministic  vector.  On  the  other  hand,  the 
clutter  in  the  range  cell  under  test  can  be  considered  as  a  large 
collection  of  point  scatterers  producing  incoherent  reflections 
of  the  radar  signal.  Then,  x  is  a  zero-mean  complex  Gaussian 
random  vector  with  covariance  matrix  E.  The  noise  e  is  a 
zero-mean  complex  Gaussian  random  vector  with  covariance 
matrix  a  1m ,  where  1m  is  the  M  x  M  identity  matrix.  In 


addition,  we  assume  that  the  clutter  reflections  and  the  thermal 
noise  are  statistically  independent. 

The  radar  dwell  consists  of  a  series  of  pulses  that  can  be 
seen  as  “snapshots”  of  the  range  cell  under  test.  If  the  pulse 
duration  is  short  with  respect  to  the  dynamic  of  the  target  and 
its  environment,  it  is  reasonable  to  assume  that  their  scattering 
coefficients  are  constant  during  each  pulse.  However,  from 
pulse  to  pulse,  we  consider  the  clutter  scattering  coefficients 
as  independent  realizations  of  the  same  random  process.  Then, 
the  distribution  of  each  snapshot  is 

yd  ~  CJ\f(AfjL,  AY>Ah  +  cr/M),  d  =  1, . . . ,  D,  (10) 

where  D  is  the  total  number  of  snapshots  or  pulses  in  the 
radar  dwell  and  H  denotes  the  conjugate  transpose  operator. 


C.  Known  and  Unknown  Model  Parameters 

We  assume  the  system  response  matrix  A  is  known,  since 
we  consider  that  the  receiver  antenna  array  has  been  calibrated. 
In  addition,  we  assume  that  the  power  of  the  thermal  noise 
a  is  known,  because  it  can  be  easily  estimated  from  the 
recorded  data  when  no  signal  has  been  transmitted.  However, 
we  suppose  that  we  have  no  prior  knowledge  about  the  target 
and  the  clutter,  nor  do  we  count  on  a  secondary  data  set  for 
estimating  the  statistical  properties  of  the  clutter.  Hence,  the 
vector  /x  and  the  matrix  E  are  the  unknown  parameters  of  the 
statistical  data  model  (10). 


III.  Detection  Test 


We  aim  to  decide  whether  a  target  is  present  or  not  in  the 
range  cell  under  test,  based  on  the  recorded  data.  Then,  the 
decision  problem  consists  of  choosing  between  two  possible 
hypotheses:  the  null  hypothesis  Ho  (target-free  hypothesis)  or 
the  alternative  hypothesis  Hi  (target-present  hypothesis) 


(  Ho  :  =  0,  E 

\  Hi  :  /i  7^  0,  E 


(11) 


where  the  matrix  E  is  considered  as  a  nuisance  parameter. 
Next,  we  derive  the  generalized  likelihood  ratio  (GLR)  test 
[7]  and  study  its  performance.  We  omit  derivation  details  in 
this  paper  due  to  space  limitation;  refer  to  [10]  for  an  extended 
version  of  this  paper. 


A.  GLR  Test 

The  logarithmic  GLR  test  decides  Hi  if 

InLGLR  =  \nf1(y1,...,yD-y1,t1)- 

ln/o(yi)  •  •  •  ,yD'>  ^o)  >  7)  (12) 

where  /o  and  /i  are  the  likelihood  functions  under  Ho  and 
Hi,  Eo  and  Ei  are  the  MLEs  of  E  under  Ho  and  Hi,  jl1 
is  the  MLE  of  /jl  under  Hi,  and  7  is  the  detection  threshold. 
For  simplicity  of  notation  we  will  omit  references  to  the  data 
in  the  arguments  of  the  functions  fo  and  fi  in  the  rest  of  the 
paper. 

We  start  the  derivation  of  the  GLR  test  by  determining  the 
likelihood  functions  and  MLEs  of  the  unknown  parameters. 
Under  hypothesis  Ho ,  it  is  assumed  that  /x  =  0;  then 

ln/0(£)  =  —D  [Mlrnr  +  In  \C\  +  tr(C_1S'0)]  ,  (13) 


where  C  =  AY>AH  +  gIm  is  the  theoretical  covariance  matrix 
of  the  data,  defined  in  (10),  and  So  is  the  sample  covariance 
matrix 

1  D 

S0  =  pJ2ydyd-  (14) 

U  d=  1 

The  MLE  of  E  is 


So  =  A+S0A+H  -  a(AHA)~\  (15) 


where  A+  =  (AH A)  1 AH  is  the  generalized  matrix  inverse. 
Under  hypothesis  Hi,  the  likelihood  function  is 


\nfi(fi,T,)  =  —D  Mlnir  +  In  |Cj  +  tr(C  1C\) 


(16) 


where 

1  D 

Cl  =  J2  (yd  -  AP>  (Vd  -  Av)H  ■  (17) 

U  d=  1 

The  MLE  of  the  unknown  parameters  are 

Ai  =  A+y ,  (18) 

Si  =  .r.S',,!"  —  a(AHA)  \  (19) 


where  y  is  the  sample  mean  vector 

1  D 

y  =  pJ2vd’  (2°) 

U  d=  1 

and  Si  is  the  sample  covariance  matrix 
1  D 

Si  =  p  (yd  -  y)  -  v)H  •  (21) 

U  d=  1 


Then,  substituting  the  MLEs  and  likelihood  functions  in  (12), 
and  after  some  mathematical  simplification,  the  GLR  test  can 
be  written  as 


Tglr  = 


1  +  yH  A  (AhS1A)  1  Ah »-• 


1  d 


(22) 


Since  (22)  is  a  monotonically  increasing  function  of  the  second 
term  inside  the  brackets,  an  equivalent  detection  test  statistic 
can  be  defined  as 

Tglr  =  VHA  (AhS1A) _1  AHy.  (23) 


B.  Detection  Performance 

Applying  Corollary  5.2.1  from  [8],  it  is  straightforward  to 
verify  that  the  detection  statistic  (23)  is  distributed  as  follows 
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Fr2P,2(D-P)  (A)  under  Hi 


(24) 


where  PUl,u2  denotes  an  F  distribution  with  vi  and  degrees 
of  freedom,  and  T' Vl,v2  (A)  denotes  a  non-central  F  distribu¬ 
tion  with  vi  and  degrees  of  freedom  and  non-centrality 
parameter  A.  The  non-centrality  parameter  is  given  by 


A  =  2  Dfi 
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(25) 


Thus,  the  detection  performance  becomes 


PFA  Q J72P,,2{D-P)  (T  ) 

-Pd  =  Qf'2p,2(d-p)W  (70)  (26) 


where  Q  is  the  right- tail  probability  function  [7,  Chap.  2]  and 
y  is  the  detection  threshold  for  the  required  probability  of 
false  alarm.  In  particular,  note  that  the  expression  for  Pp a  does 
not  depend  on  the  covariance  of  clutter  and  thermal  noise,  nor 
on  the  transmitted  signal;  hence,  equation  (23)  is  a  CFAR  test. 

C.  Detector  Performance  in  Inhomogeneous  Clutter 

A  significant  characteristic  of  the  detector  presented  in  this 
paper  is  that  the  test  statistic  (23)  is  computed  using  only  data 
from  the  cell  under  test.  No  assumption  is  required  about  the 
range  homogeneity  of  the  clutter.  Thus,  the  detector  is  robust 
against  spatially  fluctuating  clutter. 

To  support  this  assertion,  the  Pp a  is  computed  using  5  •  105 
runs  of  Monte  Carlo  simulations  in  the  presence  of  inhomo¬ 
geneous  clutter.  The  clutter  echoes  are  assumed  to  follow  a 
compound-Gaussian  distribution  with  covariance  matrix  r£. 
The  texture  r  has  a  generalized  Gamma  probability  density 
function  [4]: 

r(r'l  =  wj(d"T^'eKp(-r)'  (27) 

where  v  is  the  order  parameter  and  S  is  the  average  power. 
Note  that  v  =  00  corresponds  to  clutter  echoes  with  Gaussian 
distribution,  but  when  v  decreases,  they  deviate  from  Gaus¬ 
sian.  Additional  information  about  the  simulation  setup  is  as 
follows.  The  transmitted  signals  are  rectangular  pulses,  which 
are  transmitted  in  a  sequence  of  alternating  vertical  (V)  and 
horizontal  (H)  polarization  (K  =  2).  The  radar  dwell  consists 
of  10  repetitions  of  this  polarimetric  sequence  (D  =  10). 
On  the  receiver  side,  samples  of  the  V  and  H  electric  field 
are  recorded  ( Q  =  2)  at  a  rate  of  one  sample  per  pulse 
(N  =  1).  Since  the  simulations  are  intended  to  compute  the 
PpAi  no  target  is  considered.  The  average  power  of  the  texture 
is  S  =  50  and  the  covariance  of  the  speckle  component  is 

S  =  /3. 

For  comparison  purposes,  we  also  compute  the  performance 
of  two  other  polarimetric  detectors  that  use  secondary  data: 

•  Polarization- space-time  GLR  (PST-GLR):  This  detector 
was  derived  under  the  assumption  of  homogeneous  Gaus¬ 
sian  clutter  [1]. 

•  Texture  free  GLR  (TF-GLR):  This  detector  was 
formulated  assuming  that  the  clutter  follows  a 
compound-Gaussian  distribution  [4]. 

In  Fig.  1,  the  Pfa  of  our  detection  test  Tglr,  as  well  as  of 
the  PST-GLR  and  TF-GLR  tests,  is  plotted  as  a  function  of  the 
order  parameter  of  the  clutter  texture  distribution.  For  the  latter 
two  tests,  we  consider  80  adjacent  range  cells  for  generating 
the  secondary  data.  We  observe  that  the  performance  of  our 
detector  Tglr  and  the  TF-GLR  test  remains  constant  at  the 
designed  false-alarm  rate.  However,  there  is  a  discrepancy 
between  the  design  and  the  actual  clutter  that  increases  the 
probability  of  false  alarm  for  the  PST-GLR  test  as  the  clutter 
process  departs  from  being  Gaussian. 

D.  Detector  Performance  Using  Real  Data 

To  evaluate  the  performance  of  the  polarimetric  detectors 
in  a  real  application  situation,  we  use  data  collected  at  the 


Fig.  1.  Probability  of  false  alarm  as  a  function  of  the  order  parameter  v  of 
the  clutter  texture  distribution. 


HH  Channel 


Time 


Fig.  2.  Magnitude  in  range  and  time  domain  of  the  IPIX  radar  dataset 
stare  1,  collected  on  November  11,  1993.  The  target  location  is  indicated  by 
the  markers  “>”  and 


Osborne  Head  Gunnery  Range  (OHGR),  Dartmouth,  Nova 
Scotia,  Canada,  with  the  McMaster  University  IPIX  radar  [9] . 
Specifically,  we  use  the  data  recorded  on  November  11,  1993. 
Fig.  2  shows  the  magnitude  of  the  echoes  from  a  small  target  in 
inhomogeneous  sea  clutter  of  the  four  polarimetric  channels: 
HH,  VH,  VV,  and  HV. 

Fig.  3  shows  detection  maps  for  the  Tqlr,  PST-GLR  and 
TF-GLR  tests,  setting  the  threshold  for  Pfa  =  10-3.  For  the 
latter  two  tests,  a  guard  region  of  two  range  cells  on  either 
side  of  the  cell  under  test  is  used  to  avoid  the  leakage  of 
the  target  energy  into  the  secondary  data.  We  observe  that 
our  detection  test  finds  the  true  position  of  the  target  with 
only  a  few  false  detections  despite  the  strong  clutter  echoes, 


Fig.  3.  Detection  map  in  range  and  time  domain  for  the  IPIX  radar 
dataset  stare  1,  collected  on  November  11,  1993.  Black  pixels  indicate  that  the 
detection  statistic  is  larger  than  the  threshold.  The  target  location  is  indicated 
by  the  markers  “>”  and  “<”. 


which  can  be  appreciated  from  Fig.  2.  On  the  other  hand, 
the  PST-GLR  has  high  false-alarm  rate  due  to  the  presence 
of  inhomogeneous  clutter,  as  was  predicted  in  the  previous 
subsection.  The  TF-GLR  test  greatly  reduces  the  false-alarm 
rate  with  respect  to  the  PST-GLR  test;  however,  it  still  remains 
considerably  high. 

We  acknowledge  that  comparing  the  TF-GLR  test  with  the 
other  two  is  not  strictly  fair,  since  the  closed  form  of  this  test 
was  developed  for  only  two  polarimetric  channels  (we  have 
used  HH  and  VV  for  the  previous  results).  Our  goal  is  not  to 
extend  the  work  in  [4]  but  to  show  that  polarimetric  detection 
tests  using  secondary  data  may  have  larger  false  false-alarm 
than  their  designed  value,  in  the  case  of  inhomogeneous 
clutter. 


IV.  Target  Detection  Optimization 


We  aim  at  improving  target  detection  by  optimizing  the 
design  of  our  system.  In  the  previous  section,  we  have  shown 
that  the  target  probability  of  detection  depends  on  the  system 
characteristics  through  the  non-centrality  parameter  A,  which 
in  turn  depends  on  the  system  response  A.  We  recall  that 
matrix  A  carries  the  information  of  the  transmitted  waveform 
and  the  receiver  sensor  array.  Our  optimization  approach 
consists  of  designing  the  matrix  A  in  order  to  maximize  the 
parameter  A,  and  consequently  the  probability  of  detection. 

In  order  to  find  the  value  of  the  matrix  A  that  maximizes 
the  parameter  A,  we  rewrite  Equation  (25)  as 
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Maximizing  A,  given  /x  and  X,  is  equivalent  to  minimizing 
the  second  term  of  (28).  Assuming  that  rj  is  a  vector  whose 
entries  are  the  parameters  of  the  transmit  waveforms,  the 
system  response  matrix  is  parameterized  as  A  =  A{rj).  Thus, 


to  improve  detection,  we  seek 
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We  first  solve  the  product  AH  A  using  the  system  response 
matrix  given  by  (9),  obtaining 

K 

AHA  =  cPsYJZkh,  (30) 

k= 1 

where  the  value  of  the  coefficient  c  depends  of  the  receiver 
sensor  array,  and  Vs  is  the  transmit  signal  energy: 

N 

VS  =  J2  Ht)\2.  (31) 

t=  1 

The  constant  c  for  a  conventional  V/H  sensor  array  is  c  =  1; 
and  for  the  electromagnetic  vector  sensor  is  c  =  2.  Thus, 
the  performance  of  our  polarimetric  detector  does  not  depend 
on  the  pulse  envelope  of  the  waveform  but  rather  on  the 
polarization  parameters  rj  =  [aq,  /?i, . . . ,  olk ,  Pk\- 

We  mention  here  that  in  a  real  application  the  true  values 
of  /x  and  E  are  not  known.  Instead,  their  estimates  fi1  and  Ei 
should  be  used.  Nevertheless,  solving  Equation  (29)  with  the 
true  target  and  clutter  values  provides  an  upper  bound  of  the 
detection  improvement. 


A.  Simulation  Example 

In  order  to  study  the  detection  improvement  by  optimally 
selecting  the  signal  polarization,  we  design  a  simple  problem. 
We  assume  a  system  in  which  the  first  pulse  is  vertically 
polarized  and  the  second  can  be  of  any  arbitrary  polarization. 
When  the  second  pulse  is  horizontal,  this  system  corresponds 
to  a  conventional  polarimetric  radar.  Fig.  4  shows  the  average 
detection  performance  of  1000  cases  where  parameters  /x  and 
E  were  randomly  selected.  The  figure  shows  the  average 
performance  of  the  conventional  system  (the  second  pulse 
has  fixed  H  polarization)  in  addition  to  the  results  obtained 
when  the  polarization  is  optimally  selected  to  satisfy  Equation 
(29)  and  when  the  worst  possible  polarization  is  selected.  We 
observe  that  the  performance  of  the  conventional  system  is 
somewhere  in  between  the  optimal  and  worst  polarization. 
Specifically,  the  optimal  polarization  shows,  for  Pd  =  0.5, 
a  reduction  of  1.5dB  of  target-to-clutter  ratio  with  respect  to 
the  conventional  system  and  3.8dB  with  respect  to  the  worst 
polarization. 


V.  Conclusions 

In  this  paper,  we  addressed  the  problem  of  detecting  small 
and  slowly  moving  targets  in  heavy  inhomogeneous  clutter 
by  exploiting  polarization  diversity.  To  improve  the  detection 
performance,  we  first  developed  a  new  polarimetric  detector 
that  is  robust  against  space  and  time  variations  of  the  clutter. 
The  proposed  detector  decides  the  presence  of  a  target  based 
only  on  the  data  recorded  from  the  cell  under  test,  without 
requiring  secondary  data  or  prior  knowledge  about  the  clutter. 
We  tested  this  detector  with  synthetic  and  real  data,  and  it 
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Fig.  4.  Detection  performance  of  the  GLR  test  statistic  Tqlr  as  a  function 
of  the  target-to-clutter  ratio,  for  PpA  =  10-3. 


demonstrated  a  significant  performance  improvement  with  re¬ 
spect  to  other  polarimetric  detectors.  An  additional  advantage 
of  our  detection  algorithm  is  that  it  is  computationally  less 
intense  because  it  does  not  require  recording  and  processing 
of  training  data.  Furthermore,  we  proposed  an  adaptation  rule 
for  designing  the  transmit  polarization  by  maximizing  the 
non-centrality  parameter  of  the  detector  distribution  in  order 
to  improve  the  target-detection  performance.  We  showed  that 
the  detection  performance  depends  on  the  signal  polarization, 
but  no  on  the  pulse  envelope.  A  numerical  experiment  was 
presented  to  demonstrate  the  performance  improvement  that 
could  be  obtained  from  the  adaptation  rule. 
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